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基于 热 啊 应 法 的 航天 颖 推进 入 质量 测量 热 模型 
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摘 要 : 本 文 以 热 响应 法 测量 航天 器 微 重力 条 件 下 贮 箱 推进 剂 剩余 质 量 为 背景 ， 建 立 了 航天 器 贮 箱 内 外 热 环境 耦合 作用 下 的 整 
体 热 分 析 模型 ， 通 过 将 航天 器 贮 箱 外 部 热 环 境 视 为 第 二 类 浮动 热 边 界 条 件 ， 实 现 贮 箱 气 液 两 相 分 布下 的 热 分 析 解 耦 计算 ， 为 热 
响应 法 提供 精确 的 温度 场 计算 方 法 。 采 用 该 方法 ， 针 对 热 响 应 法 测量 微 重力 条 件 下 某 航 天 器 贮 箱 内 部 推进 剂 质量 所 需 的 温度 场 
分 布 ， 通 过 数值 仿真 获得 了 空间 在 轨 阶 段 ， 热 响应 法 加 热 工 作 时 贮 箱 内 外 热 环 境 整体 耦合 下 的 温度 场 分 布 ， 并 依据 特定 检测 点 
的 瞬 态 温度 变化 ， 反 演 得 到 了 剩余 推进 剂 的 质量 。 研 究 发 现 采 用 热 响应 法 测量 推进 剂 质量 时 ， 贮 箱 温 度 场 不 仅 受 贮 箱 内 部 加 热 
影响 ， 在 轨 外 部 热 环 境 也 会 明显 影响 贮 箱 壁 面 温 度 的 均匀 性 。 
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A Thermal Model for Measurements of Aircraft Propellant Mass Based on the 


Thermal Propellant Gauging System Method (TPGS) 


HU Zhen-wen!, SUN Ke-xing’, AI Qing!, XING Dong', ZHANG Gao-xiong” XIA Xin-lin! 
(1. School of Energy and Power Engineering; Harbin institute of Technology, 15000, Harbin 
2. Shanghai Institute of Satellite Engineering, 201109, Shanghai) 
Abstract: This paper is based on Thermal Propellant Gauging System method (TPGS) to measure propellant mass in microgravity. A 
thermal analysis model for aircraft propellant system coupled effects of space environment. Taking the effects of the outside 
environment as Neumann boundary condition. Decouple the calculations of thermal analysis and distribution in phase of gas and liquid. 
Providing a methodology for accurate thermal analysis of Thermal Propellant Gauging Method. Utilized this method analyzed a 
temperature field of propellant system. Utilized numerical simulations to calculate the heating process of Thermal Propellant Gauging 
Method and obtain the temperature distribution of propellant system. According to the temperature response of monitor point, the mass 
of propellant was calculated. It was found that the outside environment has a vital effects on thermal distribution of propellant tank. 
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0 引言 质量 已 有 成 熟 的 技术 ， 然 而 对 于 微 重力 条 件 下 ， 由 
于 其 特殊 的 环境 限制 ， 液 体 体积 测量 相当 困难 。 由 

随 着 我 国航 天 技术 的 不 断 发 展 ， 不 少 液体 推进 ” 于 液体 处 于 微 重力 或 零 重 力 条 件 下 ， 常 规 的 依赖 于 
剂 得 到 了 广泛 应 用 。 而 推进 剂 的 多 少 直接 与 航天 器 8 力 加 速度 和 浮力 的 称 重 法 ， 静 压 差 法 等 均 不 能 应 
的 寿命 相关 联 ， 故 在 航行 时 准确 的 估算 出 贮 箱 内 推 目 。 又 由 于 微 重力 条 件 下 ， 液 体 的 表面 张力 为 主导 
进 剂 的 剩余 质量 显得 至 关 重要 "1, 在 地 面 上 测量 液体 ”作用 力 ， 容 器 内 部 有 各 种 复杂 结构 ， 使 得 液体 的 位 
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置 不 稳定 ， 使 得 辐射 法 ， 光 学 法 等 地 面 测量 技术 也 
难以 应 用 外。 从 19 世纪 60 年 代 开始 , 已 经 提出 了 十 
几 种 在 轨 条 件 下 推进 剂 的 测量 的 检测 技术 PY。 其 中 ， 
热 响 应 法 在 推进 剂 剩余 量 较 低 的 时 候 上 共有 非常 高 的 
精度 ， 在 航天 器 寿命 末期 可 以 确定 推进 剂 准确 的 剩 
余 量 由。 热 响应 法 是 以 微 重力 条 件 下 推进 剂 瞬 态 热 响 
应 来 反 演 推进 剂 质 量 的 一 种 方法 ， 其 测量 机 理 主要 
是 在 测量 过 程 中 向 贮 箱 内 输入 特定 热流 ， 不 同 的 推 
进 剂 剩余 质量 在 贮 箱 内 的 气 液 分 布 不 同 。 因 此 在 输 
入 特定 热流 时 ， 其 热 响应 不 同 。 根 据 这 个 差别 ， 进 
行 反 演 计算 得 到 推进 剂 的 剩余 质量 。 传 统 热 啊 应 法 
分 析 只 考虑 微 重力 下 的 气 液 分 布 对 贮 箱 热 响 应 的 影 
啊 。 实 际 在 轨 条 件 下 ， 航 天 器 和 贮 箱 不 停 的 与 外 界 
环境 进行 热 交 换 ， 复 杂 的 热 作用 势必 影响 贮 箱 表 面 
的 热 啊 应 。 

本 文 针 对 热 响应 法 的 需求 ， 建 立 了 卫星 贮 箱 内 
外 热 环 境 辜 合作 用 下 的 整体 热 分 析 模型 ， 通 过 将 卫 
星 贮 箱 外 部 热 环境 视 为 第 二 类 浮动 热 边界 条 件 ， 不 
仅 考 虑 了 贮 箱 内 的 气 液 分 布 对 热 啊 应 的 影响 ， 同 时 
考虑 了 航天 器 在 轨 条 件 下 外 界 环 境 热 作 用 对 贮 箱 的 


影响 。 采 用 该 方法 ， 针 对 热 响应 法 测量 空间 微 重力 
条 件 下 某 航 天 器 贮 箱 内 部 推进 剂 质 量 所 需 的 温度 场 
分 布 ， 通 过 数值 仿真 获得 了 空间 在 轨 阶 段 ， 热 响应 
法 加 热 工作 时 贮 箱 内 外 热 环 境 整 体 耦 合 下 的 温度 场 
分 布 。 


1 理论 模型 
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为 了 得 到 贮 箱 准确 的 热 响 应 ， 首 先 要 计算 微 重 
力 条 件 下 不 同 推进 剂 剩余 质 量 的 气 液 两 相 分 布 。 其 
次 要 建立 在 轨 条 件 下 航天 器 内 部 与 外 界 环 境 复杂 的 
热 交 换 模型 。 最 后 将 在 轨 条 件 外 界 换 热 作为 边界 热 
流 与 微 重力 条 件 下 的 气 液 分 布 耦合 计算 ， 模 拟 热 响 
应 法 加 热 阶 段 的 热 响应 。 耦 合计 算 流 程 图 如 图 2 所 


示 。 
贮 箱 整体 耦合 热 分 析 
贮 箱 外 部 热 环境 贮 箱 内 部 热 环境 
六 热 | | 推进 剂 空间 分 布 
太 | | 球 || 自 | | 发 | | 公 | | 克 | | 由 
阳 | | 反 || 身 | | 热 -| 全 | | 体 | | 学 l 
看 | | 照 | | 红 | | 设 | | 没 | | 等 | | 扒 
照 | | 辐 | | 外 | | 备 | | 恤 | | 结 | | 型 进 
射 | | 加 0 ag 
pg 速 | | 剩 
度 | | 祭 
_ 质 
外 浮动 热 边界 年 
邮箱 整体 温度 分 布 


2 耦合 计算 流程 图 


Fig. 2 The diagram of coupled computation process 


在 轨 条 件 下 ， 航 天 器 内 外 部 件 与 外 太空 环境 进 
行 复杂 的 热 交 换 ， 传 热 模 型 如 图 1 所 示 。 
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图 1 在 轨 条 件 下 的 换 热 模 型 


Fig. 1 The schematic diagram of heat transfer model when 


in the orbital motion 


两 相 流动 的 基本 控制 方程 是 动量 方程 四: 
FPWrv (p= vp+v | pv rog+F (1) 


采用 FLUENT 进行 模拟 计算 微 重力 条 件 下 推进 
剂 气 液 两 相 分 布 。 其 中 采用 的 是 有 限 体 积 法 中 的 
VOF 计算 模型 。 过 引入 q 相 的 体积 分 数 变量 w, ， 计 
算 相 界面 。 追 踪 相 界面 主要 是 依靠 求解 单 相 或 多 相 
的 容积 比率 的 连续 续 方 程 来 完成 的 。 对 第 q 相 ， 方程 
如 下 9; 


1| 5 < 
S| le 0,piv)= 和 + 向 ) (2) 
式 中 ， 志 ,为 q 相传 递 到 p 相 的 质量 ，m ,为 p 相传 
递 到 q 相 的 质量 ， 源 项 % 为 0。 主 相 容积 比率 的 计 
算 基 于 如 下 的 约束 : 

> (0,)=1 G) 


g=1 


而 在 轨 条 件 下 航天 器 与 外 界 环境 的 换 热 计算 采用 


SINDA/FLUENT 进行 模拟 计算 。 通 过 设 定 在 轨 工 况 
和 轨道 参数 ， 建 立 航 天 器 在 轨 条 件 下 的 热 模型 ， 计 
算 热 量 分 布 。 采 用 的 计算 模型 主要 为 热 网 络 法 和 蒙 
特 卡 洛 法 。 热 网 络 法 基本 控制 方程 为 节点 的 能 量 守 
恒 方 程 : 


peF =AVT+0, +0Q,+0, (4) 


式 中 ，2. 为 表面 间 的 对 流 换 热量 ，0, 为 表面 间 的 辐 
射 换 热量 ， 0, 为 发 热 设 备 的 内 热源 项 。 蒙 特 卡 洛 法 
计算 大 型 复杂 系统 换 热量 方程 为 : 


D, ;=0(4eT :RD, , 


-AjejTy :RD,,) (5) 


式 中 , ;是 和 44 向 4; 的 辐射 传 热 热流 量 , [W] ; 4、 


' 
4 是 面积 ，[m” ] ;RD; ; 、RD, ;分 别 是 怒 对 4)、 
有 4; 对 4 的 辐射 传递 系数 

得 到 外 界 对 贮 箱 的 热 影响 以 边界 热流 的 形式 输 
出 。 之 后 的 耦合 计算 依然 采用 的 是 ANSYS FLUENT 
软件 进行 模拟 计算 ， 得 到 贮 箱 内 特定 点 的 温度 响应 。 


2 微 重力 下 推进 剂 气 液 分 布 数值 模拟 


在 数值 模拟 计算 时 ， 首 先 要 建立 贮 箱 内 部 的 几 
何 模型 和 网 格 文件 。 图 3 为 贮 箱 系统 的 几何 模型 。 
微 重力 条 件 下 液体 的 表面 张力 为 主导 ， 为 了 保证 推 
进 剂 的 持续 输出 ， 贮 箱 内 一 般 会 采用 推进 剂 管理 装 
置 进行 管理 。 推 进 剂 管理 装置 对 贮 箱 的 气 液 分 布 影 
响 显 著 ， 故 贮 箱 模 型 内 也 包含 了 简单 的 管理 装置 。 
采用 ANSYS CFD ICEM 建立 非 结 构 化 网 格 文件 , 如 
图 4 所 示 。 为 保证 计算 结果 进度 ， 网 格 密度 较 高 ， 
网 格 体 单元 数 为 818394， 网 格 面 单元 数 1689015， 
网 格 节点 数 为 164340。 
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3 贮 箱 系 统 几 何 模型 


Fig. 3 The geometry model of Propellant Gauging System 


4 贮 箱 网 格 划 分 图 
Fig. 4 The mesh of the tank 
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图 5 推进 剂 剩余 质量 60kg 时 贮 箱 内 气 液 两 相 分 布 
Fig. $ The distribution of liquid-phase volume fraction in the 


tank when the remaining mass of propellant is 60 kg 
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图 6 推进 剂 剩 余 质 量 40kg 时 贮 箱 内 气 液 两 相 分 布 


Fig. 6 The distribution of liquid-phase volume fraction in the 


tank when the remaining mass of propellant is 40 kg 


以 常用 的 液体 推进 剂 甲 基 有 愤 为 工 质 ， 模 拟 了 零 
重力 加 速度 条 件 下 , 推进 剂 剩余 质量 为 60kg 和 40kg 
时 的 两 相 分 布 。 图 5、 图 6 分 别 为 推进 剂 剩余 质量 为 
60kg 和 40kg 时 贮 箱 内 部 液 相 体积 分 数 分 布 图 , 其 中 
数值 为 1 表示 为 液 相 推进 剂 。 数 值 为 0 表示 气相 。 


日 


i pA 人 4M 不甘 日 工 || 
IV 合 作 期 证 | 


站 hinaV 
UL mInaA 


甚至 烧 坏 壁面 。 可 以 看 出 推进 剂 质量 仅 相差 20kg， 
而 检测 点 的 温度 响应 差异 显著 。 而 热 响应 法 便 是 建 
立 不 同 推进 剂 剩余 质量 检测 点 的 温度 响应 数据 库 来 
推算 实际 推进 剂 剩余 质量 ， 温 度 响 应 差距 越 大 ， 表 
明 该 方法 的 精度 越 高 ， 表 明 热 响应 法 在 推进 剂 剩余 


可 见 在 推进 剂 管理 装置 的 作用 下 ， 液 态 推进 剂 主要 
集中 在 贮 箱 下 部 份 ， 且 液 面 出 现 凹陷 ， 推 进 剂 沿 着 
导 液 模 相 上 扑 升 ， 推 进 剂 浸没 下 半球 表面 。 


3 热 啊 应 耦合 模拟 计算 


由 于 贮 箱 内 的 热 响 应 过 程 是 在 轨 条 件 下 实施 的 
热 动 作 ， 在 贮 箱 热 仿 真 过 程 中 需要 实时 考虑 航天 器 
的 影响 。 航 天 器 受到 外 热流 、 内 功 耗 的 影响 ， 贮 箱 
周围 的 环境 实时 变化 。 以 地 球 同步 卫星 为 例 ， 采 用 
SINDA/FLUENT 模拟 计算 了 在 轨 条 件 下 贮 箱 内 部 热 
环境 。 采 用 的 轨道 为 地 球 同 步 轨道 ， 计 算 夏 至 日 贮 
箱 与 环境 的 换 热 ， 得 到 在 轨 条 件 下 的 贮 箱 表面 热流 
分 布 ， 将 计算 得 到 的 热流 分 布 作为 浮动 的 热流 边界 
进行 输出 。 

提取 在 轨 条 件 下 得 到 的 贮 箱 表 面 热流 分 布 ， 与 
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之 前 计算 得 到 的 贮 箱 内 部 推进 剂 分 布 耦合 计算 ， 以 
贮 箱 表面 的 热流 分 布 为 外 界 边界 条 件 ， 计 算 在 轨 条 
件 外 界 环境 对 贮 箱 的 热 影响 。 同 时 模拟 实施 热 响 应 
的 加 热 过 程 ， 在 贮 箱 表面 特定 位 置 ， 输 入 热 响应 加 
热 热流 ， 同 时 检测 特定 检测 点 的 温度 响应 。 通 常 实 
施 热 响应 法 加 热 过程 为 10-100h 不 等 , 本 文采 用 的 加 
热 时 间 为 23n。 如 图 7 所 示 ， 为 推进 剂 剩余 质量 分 别 
为 40kg 和 60kg 时 热 响 应 法 实施 25h 贮 箱 温度 分 布 
图 。 


由 图 7 可 以 看 出 ， 在 外 界 热流 影响 下 ， 贮 箱 内 
的 温度 分 布 出 现 了 不 均匀 性 。 推 进 剂 剩余 质量 为 
0kg 时 ,温度 普遍 低 于 40kg 时 的 温度 。 通 过 检测 壁 
面 上 特定 点 的 温度 值 ， 来 反映 贮 箱 热 响应 。 图 8 为 
实施 热 响应 法 时 壁面 某 监测 点 温度 响应 曲线 。 由 图 
11 可 以 看 出 ,推进 剂 剩余 质量 为 60kg 时 检测 点 的 温 
度 均 低 于 40kg 时 的 温度 。 在 热 响应 法 整个 25h 的 加 
热 时 间 内 ， 壁 面 检测 点 的 温度 上 升 较 大 。 对 于 推进 
剂 剩余 质量 为 60kg 时 ， 最 大 温 升 为 70 度 左右 ， 而 
对 于 推进 剂 剩余 质量 为 40kg 时 ,最 大 温 升 达到 了 170 
度 。 壁 面 检测 点 的 位 置 选取 对 于 热 响应 法 十 分 关键 。 
在 推进 剂 剩 余 质量 较 少 时 ， 由 于 毛细 作用 附着 于 壁 
面 的 液体 层 厚 度 更 薄 ， 受 到 热 响 应 法 加 热 热流 影响 ， 
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质量 较 少 时 可 以 达到 很 高 的 精度 。 
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a) “推进 剂 剩余 质量 40kg 
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b) “推进 剂 剩余 质量 60kg 
图 7 热 响 应 法 加 热 25h 不 同 推进 剂 剩余 质量 贮 箱 温度 分 布 
Fig. 7 The temperature distribution in the tank when applied 


TPGS method to heat the system for 25 h. a) the remaining mass 


of propellant is 40 kg, and b) the remaining mass of propellant is 


60 kg 
——40kg 
— — 60kg 
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温度 上 升 更 大 。 特 别 要 注意 的 是 ， 加 热 片 的 位 置 要 
选 在 液体 覆盖 区 ， 若 在 气相 区 ， 温 度 会 急剧 上 升 ， 


8 热 响应 法 加 热 25h 壁面 不 同 推 进 剂 剩 余 质 量 监测 点 温度 
响应 曲线 


Fig. 8 The temperature of monitor point for different remaining 


mass of propellant when applied TPGS to heat the tank for 25 h. 


5 结论 


本 文 针 对 热 响应 法 测量 微 重力 条 件 下 贮 箱 内 推 
进 剂 剩余 质 量 ， 建 立 了 微 重力 条 件 下 同时 考虑 贮 箱 
内 气 液 两 相 分 布 和 在 轨 条 件 热 边界 的 一 体 化 分 析 方 
Ws NM SR I TE 

箱 内 部 气 液 两 相 分 布 ， 在 轨 条 件 下 外 界 环境 对 贮 箱 
的 热 影响， 以 边界 热流 的 形式 与 贮 箱 内 两 相 分 布 看 
合计 算 ， 模 拟 计算 了 热 响应 法 的 加 热 阶 段 ， 得 到 不 
同 推进 剂 剩余 质量 ， 特 定 检 测 点 的 温度 响应 。 结 果 
表明 在 轨 条 件 外 界 环境 对 贮 箱 表面 温度 分 布 有 一 定 
的 影响 ， 在 外 界 热 流 作 用 下 ， 贮 箱 表 面 温 度 分 布 出 
现 了 不 均匀 性 。 得 到 了 推进 剂 剩余 质量 分 别 为 40kg 
和 60kg 时 贮 箱 壁面 检测 点 的 温度 响应 。 推 进 剂 剩余 
质量 为 40kg 时 温度 均 高 于 60kg 时 的 温度 。 温 差 能 
够 达到 将 近 100 摄氏 度 。 推 进 剂 剩余 质量 仅 相差 
20kg， 检 测 上 温度 响应 差距 明显 。 胡 明 热 响 应 法 在 
剂 剩余 质量 较 少 时 将 有 良好 的 测量 精度 。 
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